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Abstract 

We investigate a spherical collapse model with and without the spatial curvature. We obtain 
the exact solutions of dynamical quantities such as the ratio of the scale factor to its value at the 
turnaround epoch (x = ^-) and the ratio of the overdensity radius to its value at the turnaround 
time (y = -^-) with general cosmological parameters. The exact solutions of the overdensity at the 
turnaround epoch C = Pclustcr for the different models are also obtained. Thus, we are able 

to obtain the nonlinear overdensity A = 1 + <5nl = C(|) at any epoch for the given model. The 

nonlinear overdensity at the virial epoch A v i r = Cff 2 " 1 ) i s obtained by using the virial theorem 
and the energy conservation. We obtain that the nonlinear overdensity of the Einstein de Sitter 
(EdS) Universe A^ s is 18tt 2 (^- + §) ~ 147 instead of the well known value 18vr 2 ~ 178. In 
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the open Universe, perturbations are virialized earlier than in flat Universe and thus clusters are 
denser at the virial epoch. Also the critical density threshold from the linear theory (Jj-jr at the 
virialized epoch is obtained as 7^(9-^ + 6)3 ~ 1.58 instead of tJj(127t)3 ~ 1.69. This value is same 
for the close and the open Universes. We also show that there is the minimum z ta value with 
the given cosmological parameters for each model. We find that the observed quantities at high 
redshifts are less sensitive between different models. Thus, the high redshift cluster (z v ; r ~ 0.5) 
is not a good object to probe the curvature of the Universe. The low redshift cluster (zy- n ~ 0.07 
i.e. zta ~ 0.6) shows the stronger model dependence feature. However, it might be still too 
small to be distinguished. Because we obtain the analytic forms of dynamical quantities, we are 
able to estimate the abundances of both virialized and non-virialized clusters at any epoch. Also 
the temperature and luminosity functions are able to be computed at any epoch. Thus, these 
analytic forms of x, y, and £ provide the accurate tools for probing the curvature of the Universe. 
Unfortunately, the current concordance model prefers the almost flat Universe and the above 
results might be restricted by the academic interests only. However, mathematical structure of 
the physical evolution equations of the curved space is identical with that of the dark energy with 
the equation of state u>^ e = — |. Thus, we are able to extend these analytic solutions to the general 
dark energy model and they will provide the useful tools for probing the properties of dark energy. 



Background evolution equations of physical quantities in a FRW Universe with the matter are 
given by 
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H 2 =\-A = ^Pm " -2 = ^Pcr , (0.1) 



- = ~^Pm, (0.2) 

Ob O 

Pm+3(^)p m = 0, (0.3) 

where a is the scale factor, p m is the energy densities of the matter, p CI is the critical energy 
density, and k is chosen to be +1, 0, or —1 for spaces of constant positive, zero, or negative 
spatial curvature. In terms of the ratio of the matter density to the critical density J7 m , the above 
Friedmann equation (0.1) becomes 
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which is valid for all times. 



n k = n m -i, (o.4) 



Pclustcr + 3(— )/J c lustcr = , (0.6) 



We consider a spherical perturbation in the matter density, /^luster is the matter density within 
the spherical overdensity radius R. The flatness condition is not held because of the perturbation 
in the matter (The curvature is not a constant inside the overdensity patch). Thus, we have the 
another equations governing the dynamics of the spherical perturbation [1] 
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where p c \ us ter is the energy density of the clustering matter. The radius of the overdensity R 
evolves slower than the scale factor a and reaches its maximum size i?ta at the turnaround epoch 
zt a and then the system begins to collapse. 

Cosmological parameters and the curvature of the Universe can be constrained from the growth 
of large scale structure and the abundance of rich clusters of galaxies. There have been numerous 
works related to this [2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12]. Most of them reach to the similar conclu- 
sions based on the conventional approximate solutions of the background scale factor and of the 
overdensity radius. It is natural to expect that the correct values for the virial radius and the 
nonlinear overdensity obtained from the exact solutions might be different from those obtained 
from the conventional approximate solutions. In order to investigate this, we need to obtain the 
exact solutions of physical quantities. 

Now we adopt the notations in Ref. [11] to investigate the evolutions of a and R 

x = — , (0.7) 

Ota 

y = (0.8) 
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Figure 1: i t a and r ta for the different values of z ta ,. a) -ffo^ta verses for the different values of 
z ta = 0.6 , 1.2, and 2.0 (from top to bottom), b) r ta verses Sl^ for the same values of z% a as in the 
left panel. 



where at a is the scale factor of the background evolution at z ta . Then the equations (0.1) and 
(0.5) are rewritten as 
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where dr = H(x ta .)y/fl m (x ta ,)dt = H t& ^n mta dt, Q ta = 
C = pC p^ ter lata) and £ta = s(^ta) = 1 from the equation (0.7). and represent the present 
value of energy density contrasts of the matter and the curvature term, respectively. Equations 
(0.9) and (0.10) can be solved analytically and we will obtain them. 

The analytic solution of Eq. (0.9) is given by 



dx' 



j-i 



dr 
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where F is the hypergeometric function and we use the boundary condition x = when r 
(see Appendix for details). From this equation, the exact turnaround time rt a is given by 
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#taV^Wita = Fo\/^m(l +^ta) 2 tta: 



(0.12) 



where we use the fact that x ta = 1, the relations Qta = 0^-1 0- + z t&) anc ^ T = ^taV^mta*. This 
exact analytic form of the turnaround time will be used to investigate the other quantities. 



As expected, r ta (i ta ) depends on (i.e. Ok) and z ta as given in the Eq. (0.12). We 
show these properties of r ta (i ta ) in Fig. 1. In the left panel of Fig. 1, we show i ta dependence 
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Figure 2: £ verses when we choose z ta = 0.6, 1.2, and 2.0. 

on the values of for the different choice of zta models. The solid, dashed, and dotdashed 
lines (from top to bottom) correspond to zta = 0.6, 1.2, and 2.0, respectively. Eq. (0.11) is the 
evolution of the background scale factor a and we can interpret it as the age of the Universe is 
a decreasing function of f2^. Larger implies faster deceleration, which corresponds a more 
rapidly expanding universe early on. Also larger zt a means the earlier formation of the structure 
and thus gives the smaller i ta . We also show the dependence of Tt a for the values of z ta m the 
right panel of Fig. 1. Because r ta = -fftaV^mta^ta, Tta becomes larger for the larger values of 

The evolution of y given in Eq. (0.10) described as the exact analytic solution (see Appendix) 



ArcSin[^/y] - Jy(l - y) = \/(t , when r < r ta , (0.13) 



y(l - y) - ArcSinf^y] + - = y/C(r - r ta ) , when r > r ta , (0.14) 

where r and Tt a are given in Eqs. (0.11) and (0.12). £ can be obtained from this analytic solutions 
Eq. (0.13) (or equally from Eq. (0.14)) by using the fact that yta = 1 

where we use Eq. (0.12). When ^ = 1, §, §, oj =1 and thus ( = (x) 2 - This factor (x) 2 is 
the well known value of £ for the Einstein de Sitter (EdS) Universe (f2 m = 1) [1, 13]. The general 
value of £ for open or closed Universe is given by Eq. (0.15). We show the behavior of ( in Fig. 2. 



As shown in the equation (0.15), ( is inversely proportional to r t 2 a . Thus, £ decreases as Q 



increases. Because ( is the ratio of p c iuster to p m at z = z t a, it means that the smaller overdensity 
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Figure 3: t normalized to the turnaround time for the EdS Universe verses y for close, EdS, and 
open Universes (from top to bottom). 

takes longer time to turnaround and collapse. In Fig. 2, the solid, dashed, and dotdashed lines 
correspond to z ta = 0.6, 1.2, and 2.0, respectively. After we obtain the value of £, we are able to 
obtain the values of x and y at any r without any ambiguity because the analytic solutions of x, 
y, and £ given in Eqs. (0.11), (0.13), (0.14), and (0.15) are exact. 



In order to better understand the above results, it is useful to investigate the virialized times for 
the different models. From the equation (0.14), we are able to obtain the collapsing time of each 
model normalized to the turnaround time for the EdS Universe 

t 2 / /— . „. r M \ rl 3 5 



- (y/y(l - y)- ArcSin[Vy] + tt)f [- , - , - , Qta" 1 ] , (0.16) 
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where ita,EdS is the turnaround time for the EdS Universe. We show this in Fig. 3. Perturbations 
is virialized earlier in the EdS Universe than in the close Universe, and even earlier in the open 
Universe. The dashed, solid, and dotdashed lines represent t normalized to the turnaround time 
of the EdS Universe, ita,EdS of close, EdS, and open Universes, respectively. The vertical dotted 
line depicts the ratio of the virialized radius to the turnaround one y V ir = h ■ 

If only the matter virializes, then from the virial theorem and the energy conservation we have 
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), (0.17) 



[1, 4, 5, 14]. From the above equation we are able to obtain y vlv = for any model 



where C/ C i uster = — g R is the potential energy associated with the spherical mass overdensity 

Rta 

\- (0.18) 
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After replacing y v - ir = \ into Eq. (0.14), we obtain 
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(0.20) 



x viT is obtained from Eq. (0.11) with Eq. (0.19) 
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Even though we obtain the analytic expression of x v - ir in Eq. (0.20), generally this equation can 
be solved non algebraic way. Thus, it might be useful to have the approximate analytic form 
of x v ; r for the wide ranges of cosmological parameters (f2^ and z ta ) if one analyze the form by 
interpolating x v - ir values. We use the exact form in Eq. (0.20) though. 

Before we move further for the nonlinear overdensity, we investigate the common mistake for 
A v i r in the EdS Universe. Eqs. (0.11), (0.14), and (0.15) can be simplified in EdS Universe 
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1. Then the above equations become 
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ArcSinl^y] + \ = V((r - \) 



f(xl-l) 



(0.21) 
(0.22) 
(0.23) 



where we use F[-^-, ±, 1 - 0] = 1 and r ta 
the nonlinear overdensity A c in the EdS Universe is that r c = 2r ta 

time for y c = even though one uses y v - ir = \ to obtain A c . By using this assumption, x 
from the equation (0.21). Thus, one obtains 



|. The commonly used assumption to obtain 
| which is the collapsing 
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(0.24) 



However, we know the exact relation between x and y and thus we do not need to use the 
above assumption. If we insert y v - 1T = \ in the equation (0.22), then we obtain the correct 
r vir = 1 + ^ < t c . With this correct value of T v i r , we obtain the correct x v ; r = (| + by using 
the Eq. (0.21). The correct value of the nonlinear overdensity A v j r for the EdS universe is 



A^f = c( — 

V yvir , 




147. 



(0.25) 



Thus, the minimum overdensity for the flat Universe is about 147 instead of 178. We show this in 
tables 1. As zt & increases A c approaches to 147 instead of 178. Also for the closed universe A v i r 
can be smaller than this value. From the equation (0.7) we are able to obtain the virial epoch 
from the given z t a 

Zvh = - 1 (0.26) 
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1.20 
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0.09 


5.30 


1.478 


137 


0.49 


5.37 


1.481 


140 


1.03 



Table 1: C, Avirj and z v - 1T with the given values of cosmological parameters for the open, flat, 
and closed Universes. Perturbations reach turnaround and virialization earlier in the flat Universe 
than in the close Universe, and even earlier in the open Universe. Thus, clusters are denser at 
virial epoch in flat Universe than in the close Universe and even denser in the open Universe. 



The ratio of cluster to background density at the virialized epoch z^- 1T for open and closed 
Universes becomes 
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(0.27) 



where we use Eqs. (0.3), (0.6), and (0.15). We also use the fact that y vlT = \ independent of Qta- 



From the equation (0.27), we are able to investigate the linear perturbation at early epoch 
A ^5 1 + S\i n (see Appendix for details) . 

A^l + k^°l + <5iin = 1 + I {Vq\t) 1 = 1 + jj ( VC v^^) 1 • (0.28) 

2 

It is equal to the famous result for EdS Universe (-y/C = x and S7 mta = 1), 8^ = | 3 [1]. 

Again, there is a common mistake for the value of the critical density threshold <5ii n (z v ir)- In EdS 
Universe, r c = | and one obtains 5y in (z c ) from Eq. (0.28) 

M«J = x|)» 1.69. (0,29, 

However, we show that the correct virialized epoch gives t v - 1t = (| + i)r ta from Eq. (0.19). Thus, 
the correct <5ii n (£ v ir) value is 

2 

M**) = l(V~Cfs ^ + ^) J=) 3 = 1(6 + 9vr)§ c 1.58 , (0.30) 

where we use vT = The above result given in Eq. (0.30) is true for with and without 

curvature and thus it is valid for the close, flat, and open Universes. After we obtain the critical 
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density threshold 5n n (z v i r ) we are able to obtain 5n n at any epoch by using the relation 

M*) = T^r^M^vir) (0.31) 

9 v^vir ) 

where D g is the linear growth factor. There is the exact analytic form of D g for the dark energy 
with the equation of state u>de = — 3 [15]. Mathematically, this form of D g is identical with that 
of curved space and thus we can adopt this form of D g in these models. 

From the analytic forms of dynamical quantities x, y, and (, it is straightforward to estimate 
the abundances of both virialized and non-virialized clusters at any epoch. Also the temperature 
and luminosity functions are able to be computed at any epoch [16]. Thus, these analytic forms 
provide an accurate tool for probing the effect of the curvature on the clustering. As we mentioned, 
the mathematical structure of the physical evolution equations of the curved space is identical with 
that of the dark energy with the equation of state ui^e = ~~ 3- Thus, these analytic solutions give 
the guideline for the extension of them to those of the general dark energy model and they will 
provide the useful tools for probing the properties of dark energy [17]. 



A Appendix 



First, we derive the exact solution of r as a function of x given in Eq. (0.11). After we replace 
the variables Z = ^ and T = ^, the integral in Eq. (0.11) becomes [18] 



dx' 



X2 I Ta(l - ZT)"2dT = -xzF 
Jo 3 



I 1 5 _x_ 
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ta 
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where F is the hypergeometric function and we use the gamma function relation r[l + b] = bT[b]. 

We also derive the exact analytic solution of r as a function of y given in Eq. (0.10). This 
equation is solved as [19] 



dr 2 



dy' 



We use 2 J y ' f{y")dy" = (y'- 1 to obtain 



= C 2 ± T . 



(A.2) 
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V-ci + (y'~ l 

After replacing Z = ^-y and T = the LHS of Eq. (A. 3) becomes 



(A.3) 



1 3 5 ci 
2' 2' 2^ V . 



(A.4) 
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We use the relation fyF[±, f , f , f y] = \F[\, \, |, _ a y)j p[i , | ? 3 } F [i , i , 3^2] = 
y _1 ArcSin[y], and F[a, 6, 6, y] = (1 — y)~ a to obtain 
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= c 2 i + r when r < r ta , 



ArcSin 



ci 



ciy 



ciy 
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1-^1 = c 2g - r when r > r ta , 



(A.5) 
(A.6) 



where ArcSin represents arcsine function. We choose the positive (negative) sign in t to be 
consistent with the fact that R (i.e. y) increases (decreases) as r increases before (after) it 
reaches to the maximum radius R ta ( yt a = 1) at r ta . Integral constants c\, C21, and C2 g are 
obtained from the boundary conditions ^| ta = 0, y(r = 0) = 0, and y ta = 1 
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= ± V 

ta » 
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where we use the fact that ArcSinfl] = |. r ta is obtained from Eq. (0.12) 
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From Eqs. (A. 8), (A. 9), and (A. 10), we obtain the exact analytic solution of y(r) 
ArcSin [^/y] - yy(l - y) = y/(r , when r < r ta . 
\/y(l - y) - ArcSinf^/y] + | = VC(t ~ T ta) , when r > r ta . 

" l-2y 



(A.7) 

(A.8) 
(A.9) 



(A.10) 

(A.ll) 
(A.12) 



If we use the relationships —ArcSin ^/y 



+ ? = iArcTan 



y(i-y) 
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\/y(l — y) + f > then we are able to rewrite the above equations (A.ll) and (A.12) as 

-\/y(l - y) - ^ArcTan 
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VC(t ~ Tta) , when r > r ta 



(A.13) 



(A.14) 



Thus, all of these three representations are equal to each other. 



We also check the perturbation is linear at early epoch A ^> <5ii n . We use equations (A.l) 
and (A.ll) because r < r ta in this case. Thus when we choose the limit r — ► 0, we obtain |xa ~ t 
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and |ya + Iy2 ~ y^r. Thus we obtain 

a - 1+ ^=<(^) 3t -= 5 <(^( 1 +^)) 2 = 1+ ^ - f^r)' < A - 15 > 

If we use t — i^taV^mta^ and Hta, — §^ta^ m the matter dominated epoch, then we obtain 

2 

^lin = f (v / Cv / ^mtat^) 3 - In EdS Universe (a/C = x an d ^mta = 1), this gives the famous result 

2 2 

<^iin = f (x*^) 3 • ^ or tne c l° se an d open Universe, 5u n (z v i r ) = §(\/Cf T vir) 3 = <f)(6 + 97r)3 which 
is equal to that of the EdS one. 

We thanks K. Umetsu for useful comment. 
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